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Abstract 

We study static spherically-symmetric solutions of non-Abelian gauge theory coupled to 
Conformal Gravity. We find solutions for the self-gravitating pure Yang-Mills case as well 
as monopole-like solutions of the Higgs system. The former are localized enough to have 
finite mass and approach asymptotically the vacuum geometry of Conformal Gravity, while 
the latter do not decay fast enough to have analogous properties. 

1 Introduction 

Conformal Gravity pQ (CG) was proposed as a possible alternative to Einstein gravity ("GR"), 
which may supply the proper framework for a solution to some of the most annoying problems 
of theoretical physics like those of the cosmological constant, the dark matter and the dark 
energy. 

It is therefore essential to investigate its predictions and consequences as further as possible. 
Here we choose to study static spherically-symmetric solutions of a non-Abelian gauge system 
coupled to CG. This report may be regarded as a sequel to the previous one which dealt with 
field systems with global symmetry and perfect fluid sources [2j. 

The main ingredient of CG is the replacement of the Einstein-Hilbert action with the Weyl 
action based on the Weyl (or conformal) tensor C K x^u defined as the totally traceless part of 
the Riemann tensor (we use R K \^ l , = — d^T^ + ...): 

1 R 

CkX^v = RnXfiu — -^(dKfiRxu — dniyRxn + QXuRku ~ OX^Rku) + —{Qk^OXv ~ 9kvQX^-, (1-1) 



so the gravitational Lagrangian is 



C g = -^P^C^ (1.2) 



where a is a dimensionless parameter. The gravitational field equations are formally similar to 
Einstein equations where the source is the energy-momentum tensor and in the left-hand- 
side Bach tensor W^ v replaces the Einstein tensor: 

W,u = | T, u (1.3) 
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Bach tensor is denned by: 
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Since Bach tensor is traceless, CG can accommodate only sources with Tjf = 0. The Yang- 
Mills (YM) and Yang-Mills-Higgs (YMH) systems that will be considered here are of course 
compatible with this condition. 

The general spherically-symmetric line-element may be written in terms of a single metric 
function by exploiting the conformal symmetry pQ: 



ds 2 = B(r)dt 2 - dr 2 /B(r) - r 2 (d9 2 + sin 2 9dtp 2 ). 
The non- vanishing components of Ricci tensor and the Ricci scalar are then 



p0 pr 

-TLq — JTL r 



B" B' 
~2 r~ 

and those of Bach tensor 



R-6 — Ry> 



1--B B' R _ 2(1 - B) AB' B „ 



(1.5) 



(1.6) 



W§ 



3r 4 



+ B 1 



1 1 

+ 



B" 



3r 4 3r 2 \ B 



+ 



1 1 

+ 



B" 



3r 4 3r 2 I B 



5' 



2 5' 
r B~ 



1 B'B" 1 / B" 
37 B 2 + 12 V~B" 



1 B B 



1 B 



IB' 1 



+ l B 



3r 4 



B k 



1 ( 1 

3r 4 3r 2 



B^ 
B 



B' 
+ 1 B 



/\ 2 



2 5' 
r B" 



2B' 
r ~B 



6 B 2 r B 3 B 

1 B'B" 2. f^L 
3r B 2 + 12 1 IF 



(1.7) 



1 B'B' 1 



+ 



1 B" 



6 B 2 3r B 



1 B'B' 



+ — 



1 / B" 



/ 3r B 2 12 V B 
1 B'B'" 1 B w 1 B"" 



6 B 2 



3r B 6 B 



A useful property of these components is: 

w°-K = - B < rm "" 
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In the following we will model the source by a static spherically symmetric matter distri- 
bution, that is a (traceless) energy-momentum tensor of the form T„ = diag(T® ,Tff ,Tg). 
The "inertial mass" of the matter fields will be as usual: 



M 7 = 4vr / drr 2 T${r) 
Jo 



(1.11) 



Thanks to (jl.lOfl the gravitational field equations (|1.3f) reduce to a single very simple equation: 

(rB)"" 3a 



r 2B 

which has a similar structure to the fourth order "Poisson equation" 

V 2 V 2 u = -h 
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where h(r) is the source term. In the spherically symmetric case V 2 V 2 u = (ru)"" jv and outside 
a spherical source (or in vacuum) u(r) is given by u(r) = cq + c\r + C2/V + nr 2 . The parameters 
are related to the source (assumed to extend within r < a) by 

1 f a 1 f a 

ci = — / r 2 h(r)dr , C2 = — / r A h{r)dr (1-14) 
^ Jo ",/ 

while k is free. Note that the volume integral of the matter density (i.e. of /i(r)) turns up as 
the coefficient of the linear term in the potential rather than the 1/r one. It is related to the 
fact that in this theory the potential of a point particle is linear in accord with the behavior of 
the Green function. 

Similarly, in CG the general solution around a localized spherically-symmetric source is 

B(r) = c + c\r + c 2 /r + kt 2 ; C 2 , = 1 + 3cic 2 (1-15) 

where the additional relation between the coefficients comes from the W£ = equation which 
is of a third order. We can easily express the two parameters of the exterior solutions by: 

Q rv poo 

ci = T J drr 2 (T°(r) - T^(r))/B(r) (1.16) 
ni r°° 

C2 = \ J o drr\T«(r) - T r r (r))/B(r) (1.17) 

while k is still not fixed by the source. However, in this framework k may be considered as a 
cosmological constant such that R = 4A = —12k. We notice that taking a > corresponds to 
gravitational attraction for "normal" matter with positive energy density and positive pressure. 

In the absence of a cosmological constant, the gravitational potential is asymptotically 
linear which enables one to explain the galactic rotation curves within this context [13 Q]. 

Among all the higher order gravitational theories [U [5], CG is unique in the sense that 
it is based on an additional symmetry principle. The conformal symmetry imposes severe 
limitations on the allowed matter sources. When matter is described in terms of a Lagrangian, 
it is very much constrained, but the Abelian and non-Abelian (n generators T a ) Higgs models 
are essentially still consistent with the conformal symmetry provided the scalar field "mass 
term" is replaced with the appropriate "conformal coupling" term which introduces a non 
minimal coupling to the Ricci scalar R. The matter Lagrangian which we will use here is 
therefore 

C m = lp M $)t(L>M #) _ _L Rm 2 _ ^ _ Ip^F^, (1.18) 

where = V fJj — ieA c ^T a and F% v are the n components of the Lie algebra- valued field strength 
F^ u T a . The resulting field equations are 

D At D /i $ + A|$| 2 $ + -$ = (1.19) 
D^F a ^ = -|[$trp"$) - (Z^$)tT a $] = J av . (1.20) 



The gravitational field equations are (|1.3p with 

V = T^™') + ~ (<^V A V A |<1>| 2 - V M V,|<I>| 2 - G^|<I>| 2 ) (1.21) 

rpjjntnimal) kg^g ^ e or dinary ("minimal") energy-momentum tensor of the Higgs model and 
G^y is the Einstein tensor. 
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Now we assume spherically-symmetric fields in the simplest non-trivial case of SU(2), 
namely (T a ) bc = —ie abc and $ and take the "hedgehog" and monopole forms which may 
be written most simply in terms of the spherical unit vectors in 3-space, e r , eg and e^: 

= /( r )e« 

A^dx^ = - (a(r) - 1) (e£d0 - e a e sin 0dp) (1.22) 



The components of the energy-momentum tensor are (using the ^-equation (|1.19|) and the 
monopole parametrization above): 

= ^e s + ie s „-^ + e „i + e„ 2 -£^-|)+^(/ 2 ) / (1.23) 

f 2 / R\ 1 / AB\ 
T r r = -e s + e sv + u-e vl + e v2 -^-(R --\ - — {b'+ — ) (f 2 )' (1.24) 



where we use the following abbreviations 

l R w 2 « 2 / 2 A 4 5a' 2 (1-a 2 ) 2 n _ 

e ^ = 2 S/ ' e - = ' n = I / ' 6vl = ^ ' £ " 2 = 2eM (L26) 

and the explicit expressions for the Ricci tensor and scalar should be obtained from eq. (|1.6|) . 
The field equations for the scalar and vector fields are the following second order equations 



{r 2 Bf')' (2a 2 R 
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+ 7T / " A/ d = (1-27) 



where one should again use (jl.6j) for R, and 

(B a >)> + (!- "> _ e 2 /2fl = (L28) 



Since there is only one independent metric component, it is obvious that not all the field 
equations (|1.3p are independent. Actually there is only one independent equation and we may 
use the third order one 

w; - |t; = o. (i.29) 

However, a much simpler form is again obtained by using (jl.lOp giving therefore the following 
fourth order equation for the metric component B(r): 



(rB)"" _ a 



B 
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Actually, we can rescale the variables r and / by an arbitrary length scale £ such that we get 
the dimensionless variables x = er/l and If. The coupling constants also rescale as a/e 2 and 
A/e 2 , and in terms of these we obtain the same equations as above with just substituting e = 1 
and replacing r by x (or thinking of r as dimensionless). Since there is no natural scale in the 
system due to the conformal invariance, we may use the typical length of the scalar curvature 
say, k = l/£ 2 , or as we actually did (for convenience) nl 2 = 0.1. Solutions with different values 
of k are related to each other by a simple scaling law. 
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2 Pure Yang-Mills Solutions 



Self-gravitating pure YM solutions in GR were discovered by Bartnik and McKinnon [6] for 
asymptotically flat space-time and generalized in [7J for the case of asymptotically anti-de Sitter 
(AdS) space-time. In this section we will present the analogous of the latter solutions in CG. 
We analyze the solutions with the two possibilities of positive and negative a. Indeed, negative 
value of a is a "wrong sign" choice since it yields a repulsive linear potential of localized 
solutions. However, the attractive contribution (nr 2 ) from the negative cosmological constant 
is dominant. Therefore we do not exclude this possibility. As in the GR case, gravity can 
balance in certain circumstances the gauge fields self-repulsion to allow for a globally regular 
solutions. 

The relevant set of equations is obtained by substituting / = in (jl.28p and (|l,30p . and 
they will have the following dimensionless form in the particular case under consideration : 

x{xB)"" = -3cra' 2 , x 2 (Ba')' + (1 - a 2 )a = (2.1) 

We solve the field equations with the boundary conditions: 

B(0) = 1 , B'(0) = , B"'(0) = , B"(oo) = 2k , a(0) = 1 , o(oo) = a (2.2) 

which are necessary for regular localized solutions with finite inertial mass as well as finite mass 
parameters of the fourth order gravity, c\ and C2- The constant k is positive in order for space 
time to be asymptotically AdS. Implementing these conditions in the field equations (I2.ip leads 
after some algebraic manipulations to the following asymptotic behavior of the gauge field: 

a(x) = a + — + % + . . . (2.3) 

x ar 

For analyzing the asymptotic behavior of the metric tensor we use the following parametrization 

B(x) = kx 2 + B lX + B - (2.4) 

x 

where we further define m(x) = m(oo) — Q(x)/2x. This form is very useful in getting the 
general behavior of B{x) and obtaining the following asymptotic expression 

_,, . o an — ao — 2ko2 — 4a2an + 3a?a^ + 4a2ao — a? + 8koo — 60103^ 

B(X) = KX -\ X -\ — g ± = h 

a\ 2af 
(12a 2 — 9aia3)ag — 602a 2 ag + 3(af — 4a| + 3a\a^)aQ + 2a2a^ + 12K(3aia 2 a3 — 2a| — a\a^ 1 

3af x 

aa\ 1 aa\a 2 1 a(3aia^ + 2a 2 ) 1 



10 x 3 60 



+ . . . (2.5) 



with arbitrary parameters k, ao, 04. In fact, k and ao are fixed by the boundary conditions. 

Whenever it exists, a solution corresponding to fixed a, ao and k, has definite values of the 
parameters a\, 02, 03, 04. They can be extracted from the numerical results described below. 

Along with Ref. [7] we find that the gauge function a(x) can approach an arbitrary constant 
(ao) at infinity, leading to a continuous family of intrinsically different solutions with continu- 
ously varying magnetic charge (or flux) Q m = 1 — Oq. This contrasts the asymptotically flat case 
(in GR) where the field a(x) should approach only the values a(oo) = ±1 which correspond to 
vanishing magnetic charge. Actually, it seems that the analogous asymptotically flat solutions 
in CG do not exist. Our solutions comprise therefore a two-parameter (k, ao) family for any 
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given a. The parameter k can be set to a fixed value by an appropriate scaling of the radial 
variable. 

It may be also of some interest to recall the analogous solutions of the Abelian case which 
are known in an explicit form [8]. These solutions may carry both magnetic and electric charges 
denoted p and q: 

1 Q 

-Faudx^ A dx v = -^dt Adr+p sin 9d9 A dip (2.6) 

2 r z 

and the Weyl equations lead to solutions of the form (]1.15p for the field B(r) where the relation 
between the coefficients is modified by a contribution from the magnetic and electric charges: 

B(r) = c + cir + c 2 /r + nr 2 ; eg = 1 + 3cic 2 + ^-(p 2 + q 2 ) (2.7) 

In absence of explicit solutions, we approached the above non-linear system of equations 
numerically and found that solutions exist for generic values of the parameters a and do- 

Our results are illustrated by the three attached figures. The profiles of solutions with two 
opposite values of a are presented in Fig. [TJ By inspection of the equations, it turns out 
that the sign of a is very much apparent from the behaviour of B"(x); this is indeed what we 
observe on the figure. Furthermore, the figure shows that the magnetic field as well as m(x) and 
Q(x) are also sensitive to the sign. Note that the two magnetic components, the "transverse" 
—a'/x (which we denote b(x) for further use) and the "radial" (1 — a 2 )/x 2 are very similar, 
but not identical as shown by a closer inspection. Fig. Q] also reveals that the fields reach their 
asymptotic values rather gradually unlike the Bartnik-McKinnon solutions [6] whose structure 
clearly split in three different regions of space. 

The dependence of several parameters characterizing the solutions on the central magnetic 
field 6(0) and on the magnetic charge Q m is shown on Fig. [2]- again for two opposite a's. We 
note the strong similarity between our Fig. [2}} and the corresponding plot of Ref. [J]. The 
Qm(b(0)) curves of Fig. [2^, have a similar general structure to that of Ref. [7], but an exact 
comparison cannot be made since our solutions are purely magnetic. 




log(1+x) log(1+x) 

(a) (b) 

Figure 1: Profiles of pure YM solutions with a — ±0.5. The other parameters are ao = —1 (or Q m = 0) 
and k — 0.1. (a) The gauge components and their corresponding magnetic fields. Also added is the 
Abelian field strength, (b) The metric functions: B", m(x) and Q(x). 
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b(0) Q m 

(a) (b) 

Figure 2: Several properties of pure YM solutions for two opposite values of a: a = ±0.5. K = 0.1 (a) 
Dependence on the central magnetic field 6(0). Note that all lines cross at 6(0) = which corresponds 
to the trivial solution with ao = 1. (b) Dependence on the magnetic charge Q m . 



An additional plot (Fig. [3]), shows the dependence of the solutions on a. This figure reveals 
the non linear response to a of the solutions and exhibits in particular the asymmetry between 
positive and negative a. 

A careful inspection of these three figures reveals that they exhibit a comforting agreement 
with each other. For example, the two values of 6(0) in Fig. are equal to those obtained 
from the points where Q m = (which correspond to ao = —1) in Fig. [2^. The values of c\ and 
C2 at the same 6(0) values in Fig. [2^ are equal to those in Fig. [3] taken at a = ±0.5. 

These solutions fit nicely to the general discussion in the introduction about localized so- 
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Figure 3: a-dependence of several properties of the pure YM solutions. The other parameters are 
ao = — 1 (or Q m = 0) and k = 0.1. Note that Mi is positive for all a, but the signs of both c\ and C2 
are correlated with that of a. 
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lutions in CG. The coefficient of the linear term in the asymptotic expansion (|2.5p is verified 
to be equal to c\ calculated directly from (|1.16j) . Similarly, the coefficient C2 of the 1/x term 
can also be obtained both ways. Note also that Mj is positive for all a, but the signs of both 
ci and C2 are correlated with that of a. 

It is natural to compare the solutions obtained in this section with the Abelian solutions 
mentioned above. These solutions possess both magnetic and electric charges. 

The magnetic Abelian solution is embedded in our system as the very simple solution 
a(x) = (of Eq. (]2.ip ) which has an inverse-square behavior of the field strength. The 
electric solution may be obtained just by duality. In both cases the Weyl equations lead to the 
gravitational field given by Eq. (12.7p . In particular, the expansion (|2.4p is truncated to the 
1/x term, the field B(x) presents a horizon at some x = Xh, both the metric and the magnetic 
field are not defined at the origin. 

In contrast, the non- Abelian solutions have non-trivial a(x); the coupled equations (|2.ip 
lead to the full multipole series given above for both a(x) and B(x). Generic YM (magnetic) 
fields are regular at the origin, namely have a% 7^ (as confirmed numerically). The mass 
parameter C2, is non- vanishing and was determined numerically (the dependence on a is shown 
on Fig. [3]). The next correction (the 1/x 2 term) is also non-zero for generic solutions. 

3 Monopole-Like Solutions 

Monopoles naturally emerge as topological defects in spontaneously broken theories and their 
non observation leads to some constraints on a large number of particle physics and cosmological 
models. Similarly, the degree of physical relevance of CG depends on the results of an analogous 
study. This leads us to examine if monopoles survive in this theory and, in the case they do, 
to study their qualitative properties. 

Gravitating non-Abelian monopoles (and their black holes counterparts) were studied by 
Ortiz [9], Lee et al. [10] and Breitenlohner et al. \11\ I12j both in asymptotically flat and AdS 
spaces. 

In parallel with the previous section, we obtain here the counterpart of these solutions in 
CG. In this case we get solutions with asymptotically broken symmetry if R tends to a negative 
constant (that is i? — > —12k). In fact, this system has been studied already by Edery et al. 
[151 [16] but only for a = —1/2. We go beyond these first results, addressing the domain of 
existence of solution in the (A, a) plane and studying some physical properties of the solutions; 
this needs in particular a better understanding of the asymptotic behaviour of the solutions as 
is done below. 

We solve the field equations with the usual boundary conditions for regular and localized 
solutions: 

B(0) = 1 , B'(0) = , B"'(0) = , B"(oo) = 2k (3.1) 
/(0) = , /(oo) = / 00 , a(0) = l , a(oo) = 0. 

The new characteristic here (with respect to the GR case) is that the vacuum expectation 
value of the Higgs field is not fixed by the Lagrangian, but by the asymptotic curvature namely 
foo = \J2k/\ (assuming as usual A > 0). 

Performing a detailed asymptotic analysis of the solutions with the above boundary condi- 
tions at x — > 00 leads to several possible types of behaviour. A combination of analytical and 
numerical considerations reveals that only one of these possibilities matches with the regularity 
conditions at the origin. The resulting asymptotic form is then : 
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a(x) ~ a\/x s + ■ ■ ■ 



fix) ~ /oo + + ■ ■ ■ 



B(x) ~ B + Six 9 + ra 2 + 

We note the presence of the parameter a in the exponents of the fields B(x) and f(x), but 
not in a{x). The additional term with respect to the usual quadratic term of AdS space time 
is therefore non-polynomial. For a > 0, the exponents in B,f provide a natural upper limit 
in the domain of existence of the solutions : a/A < 3/2. A similar pattern is also present in 
the case of the GR monopole [U [101 [111 02] • Note also that the condition p > (or q < 2) 
introduces a lower limit as well, a/A > —4, which is effective for a < 0. Both conditions can be 
also viewed as conditions on A for a fixed a since small enough A may violate them. We have 
limited our numerical investigation to values of the parameters a, A such that the field obeys 
the form above with a good accuracy. For instance to a S [—0.5,0.5], A £ [0.5,2]. We expect 
however the solutions to exist for arbitrarily large values of A within the above-mentioned 
domain. Finally we notice that the limit A — >• is singular even together with a — > 0. 

The asymptotic behaviour above causes the energy density to decrease only as T® ~ l/x p 
while the radial pressure decays more rapidly. Therefore, the inertial mass (jl.lip and the 
mass parameter C2, eq. (|1.17j) clearly diverge; so does the coefficient of the linear term in the 
potential ci, eq. f)l . 16j) for a < 0. For a > the integral for c\ seems to converge since then 
p > 1, but simple integration of the field equation (jl.30p shows that it actually vanishes. As 
a consequence, the solutions, although decaying to the vacuum as x — > oo, are not localized 
enough to have a finite inertial mass. This can be seen from the asymptotic behavior of the 
solutions. 



,= i 1+1+ 3 



1 



■ ^ = 2 1 3 - V 1 " lx 



5 W 1 -! 



(3.2) 
(3.3) 
(3.4) 




Figure 4: Monopole-like solutions with two opposite values of a: a = ±0.5. The other parameters are 
A = 1 and k — 0.1. (a) the profiles of the solutions; (b) demonstration of the asymptotic behavior of 
f(x), a(x) and B(x). 
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Figure 5: A-dependence of several properties of monopole-like solutions for two opposite values of a: 
a = ±0.5. k = 0.1. Note the consistency with Fig. [4] e.g. for a < 0, f(x) increases towards its 
asymptotical value, so f\ < 0, while for a > 0, f(x) decreases towards its asymptotical value, so /i > 0. 



We analyze the solutions with the two possibilities of positive and negative a. In passing 
we can note that, in the case a = 0, gravity decouples from the gauge system. Since the 
corresponding vacuum geometry is AdS space, the scalar and vector fields then constitute a 
monopole solution in an AdS background [T3j H3]. 

We have constructed the numerical solutions and studied them for both signs of a. As 
mentioned already, negative a is a "wrong sign" choice since it yields a repulsive linear potential 
of localized solutions, but we do not exclude this possibility since a negative cosmological 
constant is induced by the spontaneous symmetry breaking and the attractive nr 2 is dominant. 
The case a = —1/2 is the one considered by Edery et al. [15} I16| . The relation with their 
parameter k is a = —k/2. (Actually, there is also a difference of a \/2 factor between their / 
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and ours.) 

The profiles of typical solutions are shown in Fig. 0] for a = ±0.5. We see a significant 
difference between the two signs of a. For negative a, the field f(x) approaches its asymptotic 
value in a monotonically increasing way (accordingly, the coefficient /i is negative). For the 
opposite case f(x) crosses the value of /oo and then approaches this value from above. The 
coefficient f\ is now positive. Similarly the behaviour of the function B"(x) is affected by 
the sign. This function starts increasing from its minimal value at the origin for a < and 
approaches monotonically its asymptotic value 2k; for a > 0, B"(x) starts decreasing from its 
maximal value and attains its asymptotic value with one oscillation which is not perceptible 
on the plot. 

Fig. Hb demonstrates the asymptotic behavior given in eq. (|3.2p - (|3.4p . The powers s, p and 
q are not shown since they are given in an explicit form. Further properties of the solutions 
are presented on Figs. [5]and[6l More specifically, these figures summarize the dependence of 
the asymptotic coefficients on A and on a. 

4 Conclusion 

We have analyzed several types of spherically symmetric solutions of the SU(2) gauge theory 
coupled to Conformal Gravity: pure YM solutions and monopole-like solutions in the YMH 
system. 

To our knowledge, the pure YM solutions have not been discussed previously in the litera- 
ture. These solutions are well localized and have a well-defined "inertial mass" as well as finite 
coefficients c\ and C2 of the "exterior solution" . Solutions exist for all values of a including 
and negative ones, and for all finite ao- Their magnetic charges are therefore continuous. These 
solutions contain also the Abelian purely magnetic solutions [8] in a singular limit. 

The monopole-like solutions exhibit a much longer range behavior due to the scalar field. 
As a result, their gravitational fields do not approach the vacuum solution (|1.15p and they do 
not posses a finite inertial mass. They exist for all values of a in the range — 4A < a < 3A/2. 
These solutions are the closest possible analogues to the self-gravitating monopoles in the GR- 
YMH system. The significant differences which still exist seem to emerge mainly from the fact 
that the mechanism for symmetry breaking relies heavily on the presence of the non-minimal 
coupling to gravity. 

The kind of equations we have solved (fourth-order) is unconventional but could be treated 
with a good accuracy by our numerical methods which in this case are indispensable. 

A natural extension of this work would be to construct the black hole counterpart of the 
solutions, i.e. solutions with the metric function B(r) presenting a regular horizon at r = r^; 
that is B(rh) = 0. In the case of monopole, it would be interesting to study the dependence of 
the extremal values of a on the horizon size and to see if the domain of existence of conformal 
monopole-black-hole fits in a pattern similar to the one of Fig. 7 of Ref. 

acknowledgment 

One of us (Y.B.) thanks the Belgian FNRS for financial support. 

References 

[1] P. Mannheim, Prog. Part. Nucl. Phys. 56, 340 (2006). 

[2] Y. Brihaye and Y. Verbin, Phys. Rev. D 80, 124048 (2009). 

[3] P. Mannheim and D. Kazanas, Astrophys. J. 342, 635 (1989). 



11 



[4] H. J. Schmidt, Int. J. Geom. Meth. Mod. Phys. 4, 209 (2007). 

[5] L. Fabbri, Higher-Order Theories of Gravitation , arXiv:0806.2610 [hep-th]. 

[6] R. Bartnik and J. McKinnon, Phys. Rev. Lett. 61, 141 (1988). 

[7] J. Bjoraker and Y. Hosotani, Phys. Rev. Lett. 84, 1853 (2000). 

[8] P. Mannheim and D. Kazanas, Phys. Rev. D 44, 417 (1991). 

[9] M. E. Ortiz, Phys. Rev. D 45, R2586 (1992). 

[10] K. M. Lee, V. P. Nair and E. J. Weinberg, Phys. Rev. D 45, 2751 (1992). 
[11] P. Breitenlohner, P. Forgacs and D. Maison, Nucl. Phys. B 383, 357 (1992). 
[12] P. Breitenlohner, P. Forgacs and D. Maison, Nucl. Phys. B 442, 126 (1995). 
[13] A. R. Lugo and F. A. Schaposnik, Phys. Lett. B 467, 43 (1999). 
[14] A. R. Lugo, E. F. Moreno and F. A. Schaposnik, Phys. Lett. B 473, 35 (2000). 
[15] A. Edery, L. Fabbri and M. B. Paranjape, Class. Quant. Grav. 23, 6409 (2006). 
[16] A. Edery, L. Fabbri and M. B. Paranjape, Can. J. Phys. 87, 251 (2009). 



12 



